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Abstract 

Aiming to establish a rigorous link between macroscopic random motion (described e.g. by 
Langevin-type theories) and microscopic dynamics, we have undertaken a kinetic-theoretical study 
of the dynamics of a classical test-particle weakly coupled to a large heat-bath in thermal equilib- 
rium. Both subsystems are subject to an external force field. From the (time-non-local) generalized 
master equation a Fokker-Planck-type equation follows as a "quasi-Markovian" approximation. 
The kinetic operator thus defined is shown to be ill-defined; in specific, it does not preserve the 
positivity of the test-particle distribution function /(x, v;t). Adopting an alternative approach, 
previously proposed for quantum open systems, is proposed to lead to a correct kinetic opera- 
tor, which yields all the expected properties. A set of explicit expressions for the diffusion and 
drift coefficients are obtained, allowing for modelling macroscopic diffusion and dynamical friction 
phenomena, in terms of an external field and intrinsic physical parameters. 
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The link between macroscopic random motion and microscopic particle dynamics has been 
a long standing fundamental problem, lying in the very heart of Non-Equilibrium Statistical 
Mechanics. In a generic manner, fluctuations due to particle interactions {collisions) are 
modelled by a Fokker-Planck-type equation (FPE), typically related to a Langevin-type 
equation of motion, which may either be derived intuitively, via physical phenomenology 
or, formally, through kinetic-theoretical arguments Q]. In the latter framework, a number 
of works in Non-Equilibrium Statistical Mechanics have been devoted to the study of the 
relaxation of a small subsystem (e.g. a test-particle, t.p.) weakly interacting with a heat 
bath. A common aim of such studies is the derivation of a kinetic equation, describing 
the evolution in time of a phase-space probability density function. This is achieved by 
using either perturbation theory (typically a BBGKY hierarchy of equations for reduced 
distribution functions (d.f.) P, Isf) or formal theories for open systems (e.g. projection- 
operator methods ^). In a generic manner, both approaches rely on a generalized master 
equation (GME), which is obtained in 2nd order in the (supposedly weak) inter-particle 
interaction; this is a non-Markovian equation, i.e. is characterized by non-locality in time. 
The kernel of the GME needs to be evaluated along particle trajectories, so the influence 
of external force fields on the microscopic laws of motion is in principle expected to modify 
the form of the collision operator. 

This paper is devoted to the rigorous derivation of a Fokker-Planck-type kinetic equation 
from microscopic dynamics, taking into account the existence of an external force field and 
interactions (possibly of long range) between particles. This work aims addresses a number 
of fundamental questions arising in the description of open systems |^ and, as a matter of 
fact, complements previous studies relying on formal projection operator methods 0,Q, 0], 
where these issues were first considered with respect to classical systems. 

a. The model . We consider a test-particle (t.p.), say S, surrounded by (and weakly 
coupled to) a homogeneous reservoir i?; by X = (x, v) = (xs(t), vs(t)) and Xr, = 
{Xj} = {xj(t), Vj(t), j = 1, 2, 3, A^} we shall denote the coordinates of the test- (S— ) 
and reservoir- (-R— ) particles, respectively. Both subsystems are subject to an external 
force field, which affects the particle dynamical trajectories {x(t),v(t)}. In the absence 
of interactions (collisions), particle trajectories correspond to the solution of the (single 
particle) dynamical problem of motion (P^{t)/dt^ = Fq. 
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The Hamiltonian of the system is: 

H = Hr + Hj: + \Hi. (1) 

Here Hr {Hy) is the Hamiltonian of the reservoir (t.p.): Hr = Ylf=iHj + Y.j<nY.n=i^jn- 
The single-particle Hamiltonian Hj (j = 1, 2, A^, S) takes into account the external field. 
Hj stands for the interaction (assumed to be weak: A <^ 1) between S and R: Hj = 
J2n=i'^^n, where Vij = V{\xi — Xj|) = 1, 2, A^, S) is a binary-interaction potential 
(possibly a long-range e.g. Coulomb-type one). The resulting equations of motion are: 

i = v, v = Fo(x,v) + AFi„t(x,v;XR;t) (2) 

The force Fq is due to the external field. The interaction force Fint(x, v; Xr; t) = 
~^Z)^(|x — Xjl) (actually the sum of interactions between S— and N /2— particles) rep- 
resents a zero-mean Gaussian random process, for a reservoir in a homogeneous equilibrium 
state p. 

Let us assume that the zeroth-order ('free') problem of motion - namely Q for A = - 
yields a known analytical solution in the form: x*^°^(t) = x + /q (it' v(t') = M(t)x + N(t)v 
[so v^^^t) = M'(i(:) X -|- N'(t) v; the prime denotes differentiation with respect to time t], i.e. 

(:::;;)=(:;::,:)(:)■■»(:) 

with the initial condition {x, v} = {x^°^(0), v*^°^(0)}, implying E(0) = I (the unit matrix). 
For a given dynamical problem, say in d dimensions {d = 1,2,3), the 2d x 2d matrix E(t) 
in © satisfies the group property: E(t)E(t') = E(t + t') \/t,t' G 3?, implying E{-t) = 
E"^(t) and a number of explicit relations among the dxd sub-matrices {M(t), N(t)}, whose 
form depends on the particular aspects of the dynamical problem into consideration. This 
working hypothesis requires that a solution of the linearized problem of motion be known. 
Specific paradigms obeying Eq. Q include: 

(a) linear oscillators, for which (in ID): M{t) = N'(t) = coscuot and uJoN(t) = —cu^'^M'^t) = 
sin uJot (here ujq is the characteristic spring vibration eigenfrequency) ; 

(b) charged particles (charge q, mass m) subject to a uniform magnetic field B = Bz along 

the 2;— axis, whose (spiral) motion obeys exactly: M = I (hence M' = 0), N'(t) is the 

/ cos Qt sin Qt \ 



rotation matrix: N'(t) 



sin Qt cos Qt 
1/ 
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and N(t) = J* dt' R(t'), where n is the 



gyroscopic (cyclotron) frequency fl = qB/m P. Il^: 

(c) the free motion (Chandrasekhar 11 1 or Rosenbluth-MacDonald-Judd limit, in the 



absence of external field(s): x(i(:) = x + tv [where \(t) = v = est.], viz. M = t^^N = I. 

b. Statistical formulation. Let p = p({X, Xr.}; t) be the total phase-space distribution 
function (d.f.), normalized to unity: / c/X p = 1. The equation of continuity in phase space 
r reads: 

9p dp d , 1 ^ , ^ 
dt -^Sxj 5vj m 

where a summation over j (= 1, 2, 3, A^, S) is understood. 

Defining appropriate 's-body' {s = 1, 2, 3, ...) reduced distribution functions (rdf ), among 

which the (1— body-) test-particle rdf: /(x, v;t) = (J, p)/j = /p^ c?Xrp (normalized to 

unity), and then integrating the total ((A^ + 1)— particle) Liouville equation (jH), one obtains 

a BBGKY hierarchy of coupled evolution equations for the rdfs. This standard procedure 

P] has here been adapted to a test-particle problem, namely by defining two types of rdfs, 

depending on whether or not the t.p. is included in the particle cluster considered ^\. 

Truncating the hierarchy, one obtains 

{dt - )) /(X; t)^X' J rfXi Lig{X, Xi; t) , 
{dt - lP - 4^)) g{X, Xi; t)^\Lj F, (Xi) /(X) , (5) 

where Lq^ {j G {S, Ir}) is the "free" Liouvillian (in the presence of the field): 

and Lj = L^i is the binary interaction operator: 

Li = -Fi„t(|x - xil) - — . (7) 

\m av nil oviJ 

Recall that the "tags" E (1) refer to the t.p. (one, any, of the reservoir particles); thus, 
/ = /(X;t), Fi(Xij^) and /2(X,Xi;t) denote the S-l-body, i?-l-body and {Ir + S)-2- 
body rdfs respectively, while g = (^(X, Xi; t) is the 'two-body' {Ir + T,) correlation function: 
^(X,Xi;t) = /2(X,Xi;t) - Fi(Xi)/(X; t). Note that the mean-field (Vlasov) term (~ A^) 
cancels, for reasons of symmetry, since the reservoir is assumed to be in a homogeneous 
equilibrium state Fi = n</)eg(vi) {n = N/V is the particle density), i.e. dFi/dt = L^^^ Fi = 
0. 
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Neglecting correlations at t = 0, Eqs. (0) lead to the Non-Markovian Generalized Master 
Equation (GME): 

dtf -Lof = X^n J^dr J t/xi rfvi Lj Uo{t) Lj 0e,(vi) / (8) 

[here / = /(x, v;t)], where the "free" Liouville operator Lq = Lq was defined in (0), Lj 
is the binary interaction Liouville operator L^i [cf. ((Tj)] and Uo{t) = Uq^\t) Uq'\t) is the 
evolution operator (propagator) defined by the formal solution of the "free" (collisionless) 
Liouville equation [i.e. (jSk) for A = 0]: f{t) = e^o''*/(0) = U^^\t) /(O) (for j G {S, 1}). 

The procedure described here formally amounts to defining the projection: P ■ = 
Or J dTpi ■ , where api and F/j denote the reservoir distribution function and phase space, 
respectively (notice that: = P ), and seeking an evolution equation for the t.p. 

rdf f = o"^^ P p. This may be achieved by defining the complementary projection, say 
Q = I — P (where I is the identity operator_), and then deriving a pair of equations de- 
scribing the evolution of Pp and Q p in time (yj. This idea has been rigorously elaborated 



in Refs. 

c. A ' quasi- Markovian' approximation - the O— operator. A widely used "markovian- 
ization" procedure consists in substituting with the zeroth-order solution, i.e. assuming 
that fit — r) ^ e~^°'^ fit) = Uoi—r) fit), and then evaluating the kernel asymptotically 
i.e. taking the time integration limit t to infinity. A kinetic operator is thus formally 
defined, here denoted as the 6— operator. Let us point out that the time-propagator 
f/(t) does not commute with /"-space gradients indeed, one rigorously obtains: 

Uo\t)^U^^\-t) = Nfit) ^ + (t) ^, for J = S, 1^ (a similar expression holds 



for ^) |9|. See that the field inevitably enters - via the dynamical matrices Nj(t), ... - the 
collision term. 

Limiting ourselves to a spatially uniform system: / = /(v; t) and substituting from 
dni) and ((Zj) into the GME (jH)), the 'markovian' approximation described above leads to a 
parabolic PDE in the form 

+ -Fo/ = — A / + —TP a/ , 9 
ot m av av av mi av 

i.e. a Fokker-Planck-tjpe diffusion equation 

df 1 df d 

at m av ovi oviOVj 



(here i,j = 1, 2, 3 = Vx,Vy,Vz)- Here, D = A is a (positive definite) diffusion matrix, which 
is given by: 



roo r r 

dr J c?xi J dvi(f)eq{vi) 
xFi„t(|x(°) - xi^l) ® FU\^('\-r) - xiW(-r)l) N'^(r) 

dTC{t,t-T)N'^{T). (11) 



Jo 

Notice the exphcit appearance of the Green-Kubo coefficients (for the interaction forces 
Fint): diffusion coefficients are thus related to the force correlation matrix Cij = 
{Fint,i{t) Fint,j(t — T))ji [= Cijij), in a stationary process], in agreement with phenomeno- 
logical stochastic theories. Defining the Fourier transform Vk of V{r), this relation takes the 
form: 

D'yj^ = ^(27r)^ r dr f dv^ 0,,(vi) / dk e^"^^^ K iV'_(r) (12) 

{d = 3, in a 3D problem), where the exponent Ar = r(t) — r(t — r) (viz. r = x^^-* — Xi'-'^^) 
can be computed by making use of Q; e.g. Ar(r) = N(r)(v('^) — Vi^^^), as in paradigms (b) 
and (c) above, where M = I P]. The vector F in the right-hand-side {rfis): Ti = = 
(1 + ^) = — ^ai, represents the dynamical friction force suffered by the particle, 
due to interactions with its environment (the reservoir). 

Summarizing so far, we have derived a Fokker-Planck kinetic equation, involving an ana- 
lytical (differential) collision operator; a set of exact expressions for the diffusion coefficients 
and the dynamical friction vector are thus obtained, to be computed for a given problem 
(taking into account the external field) either analytically or numerically; cf. e.g. in Ref. 
for case (b) above. 

In the general case of an inhomogeneous d.f. / = /(x, v; t), one obtains the kinetic 
equation 

ot ax m av av av av ax mi av 
where the above definitions still hold; G is obtained from rhs i\ll^ upon N'^ ^ N^. The 
r/ts()13p may readily be cast in the form of a Fokker-Planck operator [i.e. formally as in 
rhs ^TU^ . yet now for 6 variables (in a 3D problem): i,j = 1, 2, 6 = x,y, z, v^, Vy, Vz]. This 
6+1 variable FPE now involves a 6D friction vector: JF = (0, 0, 0; JF^^), JF^^)) and a 

6x6 diffusion matrix, in the form: ' 

iG A 



Here is a crucial remark to be made. From a fundamental point of view, the action of 
a kinetic evolution operator has to preserve the reality, the norm and the positivity of the 
d.f. / (which is defined e.g. as a normalized probability distribution function); in formal 
words, a semi-group should be defined. Furthermore, an H-theorem should be satisfied 
P], dictating monotonic (irreversible) convergence towards a state of equilibrium (where an 
entropy function S = —H attains its maximum). In fact, the diffusion matrix T)^^\ involved 
in the 0— operator, is not positive definite, as can be readily checked (see that DetP^®) = 0); 
therefore (fTUI) determines an ill-defined kinetic operator, whose action does not preserve the 
positivity of the d.f. f (and neither satisfies an H-theorem). This issue has often been 
overlooked in the past, since most studies were limited to spatially homogeneous systems 
(where the problem is not manifested, since is, in principle, positive defininite). 

d. Towards a Markovian FPE - the (^—operator. The forementioned problem has been 
mown in quantum mechanics 14,^,1131, despite its being overlooked in the classical context 

nil i^M, J' p ^ |_. 

21, In his remarkable work on quantum open systems p^, E.B. Davies proposed a new 
Markovian evolution operator, which he formally showed to define a semigroup, and was 
thus psoposed as a correct description of open systems' evolution towards the equilibrium 
state of a thermostat. "Davies' devi ce" 1^] was recently successfully applied in a (classical) 
plasma kinetic-theoretical context 0, llol|. The main stepstones and outcome of this method 
will be briefly exposed in the following. 

The alternative ($— ) operator introduced by Davies (also see in Q]), essentially amounts 

to 

= hm ^ r dt'U^'\t') . U^'\-t') (14) 

T^oo Zl J-T 

(applied to the kernel of the master equation derived above); recall that U^^\t) = exp 
is the zeroth-order (coUisionless) evolution operator (propagator). 

The analytical construction of the $ operator for a given physical problem assumes knowl- 
edge of the exact form of the dynamical matrices (defined above) , in addition to the reservoir 
equilibrium state (peqi^i) (i-e. typically a Maxwellian). The computation may be quite te- 
dious, yet quite straightforward: it is, in fact, rather simplified in the case of oscillatory 
test-particle (zeroth-order) motion, as in paradigms (a) and (b) above [unlike in free mo- 
tion (c), which is unbounded and characterized by a continuum spectrum, giving rise to 
ill-defined integrals e.g. lim^^oo J^T^'^dt]; further details, here omitted for brevity, can 



be found in 



a. 



We may briefly state the outcome of tlie calculations in two typical cases, namely 
paradigms (a) and (b) above, where the $ operator was shown to provide a well-defined 
FPE, which was proven to fulfill all the necessary requirements. In both cases, the FP-type 
kinetic equation obtained features a novel collision term, which accounts for diffusion in 
real, position space, in addition to velocity space diffusion [present, alone, in the B— FPE 
fllO|) ]. Furthermore, a modified cross- velocity-position diffusion part is obtained [replacing 
the "evil" second term in rfe ()13|) ]: finally, the dynamical friction vector JF is modified [cf. 
the last terms in the two Eqs. below, which are absent in rfe fjl3|) ] 

In the case of charged particle motion in a uniform magnetic field, the long calculation 
yields a kinetic equation in the form 



•9/ , q df 

ot ax m av 
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dvxdy 



dvydx 



+ ^[^Il(v)/] 



-2 \T^{XX)^ 
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W f 







+ [D 
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+ 



92 



dx"^ dy 
d 



■)/ 



(cf. variable definitions in Sec. I (Jap : note the cylindrical symmetry in the collision term 
(rhs), as imposed by the field (~ z). Also note the XY-space diffusion term in the second 
line. The lengthy expressions for all coefficients (involving multiple integrations, which have 
been calculated a naly tically to some extent, and then computed numerically) are exposed 
elsewhere 10, ll^ . It is verified that all of the expected mathematical requirements 
are fulfilled by this new plasma kinetic operator. Interestingly, by adopting as a working 
hypothesis that Dij and J-'i/vi are constant (approximately true for low t.p. velocities), 
Eq. (fTH|) models a multi-dimensional Ornstein-Uhlenbeck random process 0, which can be 
solved exactly jisj. 

Considering a "ID gas" of linear oscillators (in a parabolic potential field ~ uJqx'^/2), we 
have obtained: 



dt 



df 

! 

dx 



df 



i2 r 



Dvv{v)f 



+ 



92 



dvdx 



Dvx{v)f 

d_ 

dv 



d 



2 r 



dx"^ 



Dxx{v)f 





d 




:Fv{v) f 


dx 


J'xiv)f 



(16) 



All (scalar) coefficients can be explicitly computed (relying on the formulae derived above), 
yet are omitted here for brevity. The collision term thus defined (by the rhs) can be shown 
to be well-defined as a kinetic evolution operator. 
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e. Conclusion - discussion. Relying on first statistical mechanical principles, we have 
discussed the fundamental issues involved in the derivation of a Fokker-Planck kinetic equa- 
tion (FPE), for the description of test-particle dynamics in the presence of an external field. 
Exact expressions were derived for the diffusion and drift coefficients, which may be exactly 
computed for a class of physical systems. We have thus exactly recovered the anticipated 
role of the Green-Kubo (force) coefficients, relating diffusion coefficients to force correla- 
tions, in agreement with phenomenological stochastic theories. The collision mechanism 
is nevertheless an intrinsic part of the formalism here, and is not plainly represented by 
ad hoc assumptions on the nature of the process. Furthermore, the external force field ap- 
pears explicitly in both correlation functions and diffusion/drift coefficients appearing in the 
Fokker-Planck kinetic equation for the test-particle. The results presented here are qualita- 
tively reminiscent of previous formal studies 0, Q , where these issues were first discussed 
with respect to classical systems. 
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